Bound State Calculations of the Three-Dimensional Yakubovsky 
Equations with the inclusion of Three-Body Forces 

M. R. HadizadetQ and S. BayegaiJll 

Department of Physics, University of Tehran, P.O.Box 14395-547, Tehran, Iran 

(Dated: April 29, 2008) 

Abstract 

The four-body Yakubovsky equations in a Three-Dimensional approach with the inclusion of the 
three-body forces is proposed. The four-body bound state with two- and three-body interactions 
is formulated in Three-Dimensional approach for identical particles as function of vector Jacobi 
momenta, specifically the magnitudes of the momenta and the angles between them. The mod- 
ified three dimensional Yakubovsky integral equations is successfully solved with the scalar two- 
meson exchange three-body force where the Malfliet-Tjon-type two-body force is implemented. 
The three-body force effects on the energy eigenvalue and the four-body wave function, as well 
as accuracy of our numerical calculations are presented. The four-body Yakubovsky equations in a 
Three-Dimensional approach with the inclusion of the three-body forces is proposed. The four-body 
bound state with two- and three-body interactions is formulated in Three-Dimensional approach 
for identical particles as function of vector Jacobi momenta, specifically the magnitudes of the 
momenta and the angles between them. The modified three dimensional Yakubovsky integral 
equations is successfully solved with the scalar two-meson exchange three-body force where the 
Malfliet-Tjon-type two-body force is implemented. The three-body force effects on the energy 
eigenvalue and the four-body wave function, as well as accuracy of our numerical calculations are 
presented. 
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I. INTRODUCTION 

The topic of three-body forces (3BFs) is as old as nuclear physics ^ and, based on meson 
exchanges, various processes have been proposed in the past (for a review see [2I). Among 
them, the Fujita-Miyazawa force 3| with an intermediate A generated by the exchange of 
two pions is most obvious and is implemented in all modern 3BF models. However the 
nature of these 3BFs is still not completely understood. In recent years there has been new 



progress in understanding t 
perturbation theory (xPT) 



le form of nuclear forces, because of the application of chiral 
^- 8| . From this developments one can expect a more system- 
atic understanding of the form of two-body (2B) and 3B forces. However, xPT implies 
a priori unknown constants, the low-energy constants, which have to be determined from 
experimental data. The bound states of few-nucleons seem to be an ideal laboratory to de- 
termine 3BF parameters, as the binding energies are sensitive to the 3B interaction and they 
are expected to be governed by the low-energy regime of nuclear physics j^, [l^. Therefore, 
the understanding of nuclear few-body bound states is an important contribution to the 
understanding of the 3BF. To this aim and for their numerical investigations one requires 
reliable methods leading to the solutions of the non-relativistic Schrodinger equation. 

In the past several solution methods have been developed and applied to the four-body 



bound state problem by using realistic nuclear potentials, the CRCGV 

" " n 



HH 



Hi, the SV pj, the 



13|, the GFMC [14], the NCSM [15|, EIHH [l&l and the Faddeev-Yakubovsky (FY) [17|- 



271]. These calculational schemes are mostly based on a partial wave (PW) decomposition. 
Stochastic and Monte Carlo methods, however, are performed directly using position vectors 
in configuration space. One of the most viable approaches appears to be the FY method. 
The calculations based on FY are performed after a PW expansion with phenomenological 
potentials in configuration space jl^, Q, and in momentum space [3]- {2^ and recently 
with chiral potentials in momentum space {25!- [2^]. 

The FY scheme based on a PW decomposition, which includes spin and isospin degrees of 
freedom, after truncation leads to two coupled sets of a finite number of coupled equations in 
three variables for the amplitudes. In PW decomposition the number of channels that must 
be included grows very rapidly in this case, and a further complication is arisen where there 
are now six spatial dimensions rather than the three required for three-body calculations. So 
in a PW decomposition one needs a tremendous number of partial waves to find converged 
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results. In view of this very large number of interfering terms it appears natural to give up 
such an expansion and work directly with vector variables. 

On this basis recently we have extended the Three- Dimensional (3D) approach, which 
greatly simplifies the two- and three-body scattering and bound state calculations without 



using PW decomposition [28|-[38(], to the four-body bound state [39|, |40|. We have for- 
mulated the Yakubovsky equations with only 2BFs as function of vector Jacobi momenta, 
specifically the magnitudes of the momenta and the angles between them. We have obtained 
two coupled three-dimensional integral equations in six variables for the amplitudes which 
greatly simplifies the calculations without using PW decomposition. The obtained three- 
dimensional integral equations have been solved successfully for simple NN force models. In 
this paper we follow the same approach and consider the 3BFs in four-body bound state 
problem. As a simplification we neglect spin and isospin degrees of freedom and study the 
four-boson bound state problem. 

So the purpose of this work is to demonstrate that one can solve the Yakubovsky equations 
for four-body bound state without using PW decomposition and in the first attempt we have 
done it by using very simple 2B and 3B model interactions. 

In our formulation we work directly with vector variables in the Yakubovsky scheme in 
momentum space. Here the dependence on momentum vectors shows that our 3D repre- 
sentation in comparison to traditional PW representation avoids the very involved angular 
momentum algebra occurring for the permutations and especially for the 3BFs and the full 



y whereas the PW representation of the amplitudes 



solution can be reached exactly and simp 
leads to rather complicated expressions [22 1 

We believe that this work is another step forward in the development of 3D approach for 
studying the few-body systems and it is the first attempt towards the solution of the 4N 
bound state problem with the inclusion of 3NFs without performing the PW decomposition. 

This paper is organized as follows. In section [III we briefly represent the coupled 
Yakubovsky equations for four-body bound state with two- and three-body interactions. In 
section HIT] we evaluate the matrix elements of 3BFs. In section HVl we discuss our choice for 
independent variables for the unknown amplitudes in the equations and in their kernels. Sec- 
tion |V] describes details of our algorithm for solving coupled Yakubovsky three-dimensional 
integral equations. In section |Vl] we present our results for three- and four-body binding 
energies with and without model 3BFs and we provide the test of our calculation. Finally 
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we summarize in section IVIII and provide an outlook. 



II. MOMENTUM SPACE REPRESENTATION OF YAKUBOVSKY EQUATIONS 
WITH 3BFS 

The bound state of the four-body (4B) system, in the presence of 3BFs, is described by 



two coupled Yakubovsky equations 24| : 



l^i) = GotP[{l + Pum) + \M + il + Got)GoWgl\^) 

IV'2) = G'otP[(l + P34)|^l) + |^2)] (1) 

where the Yakubovsky components l-i/^i) and |'!/^2) belong to "3 + r'(123,4; 12) and "2 + 
2" (12, 34; 12) partitions of the four particles respectively. Here the free four-body propagator 
is given by Gq = {E — Hq)^^, and Hq stands for the free Hamiltonian. The operator t is the 
two-body transition matrix determined by a two-body Lippman-Schwinger equation. P, P 
and P34 are permutation operators. P = P12P23 + -P13-P23 permutes the particles in three- 
body subsystem (123) and P = P13F24 interchanges the two two-body subclusters (12) and 
(34). The quantity ^^{23, as shown in Fig. [1], defines a part of the 3BF in the cluster (123), 
which is symmetric under the exchange of particles 1 and 2 and which can be related by an 
interchange of the three particles to two other parts and VF123 ^^at sum up to the total 
3BF of particles 1, 2 and 3: Wus = W^l + wi^l + wj^l The total 4B wave function |^) is 
given as: 



\^) = {1 + P + P,,P + P)[(l + P34)|V^i) + |V^2)] (2) 

The symmetry property of under exchange of particles 1 and 2, and \1lJ2) under 

separate exchanges of particles 1,2 and 3,4 guarantee that |\E') is totally symmetric. It can 

(3) 

easily be verified that the inclusion of the 3BF component VF^gs definition of the 

first Yakubovsky component I'i/'i) does not change its symmetry property. 

In this paper we follow the notation introduced in Ref. [40] and work in a 3D momentum 
space basis. According to the two types of chains (123,4; 12) and (12,34; 12) there are two 
type of basis states. Fig. [2], which are suitable to represent the two Yakubovsky components 
1^1) and \1IJ2) in the coupled equations The representation of coupled equations (P) in 
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FIG. 1: Diagrammatic representation of the part W^23 ^ two-meson exchange 3BF. Here particle 
3 is single out by the meson-nucleon amplitude described by the blob. 
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FIG. 2: Definition of the 3 + 1 and 2 + 2 type of Jacobi coordinates. 



these basis sates will be exactly the same as displayed in Ref. 40| except that an extra term 
with Vt^i23 occurs in the first component. This is 



D^u {ui 'U2'U3|(1 + Goti2)Go\u[u2u'^) 

X(S;^'2«3|W^123|*) (3) 

where D^u = d'^Ui d^U2 d^u-^. The first matrix element can be handled as described in Ref. 



40l |. The second matrix element involves the 3BF, which has been worked out in Ref. 33 1 
in a 3D momentum space basis for three-body system. After evaluating the first matrix 
element in Eq. ([3]), the coupled three dimensional Yakubovsky integral equations can be 
rewrite explicitly as: 
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X 



E 



m 4m 3m 



d u'2 {ui\ts{e)\-U2 + U2) 



{^2 + ^^2 ^2 ^l^/'l) 



+ (^2+2^2 



2 1 2 

^2 - 3^i3 2^2 - 3^31^2) 



+ <! (^1^2^3|W^1?^|*) 



d% ^!!j^-^tf^^,. (^^>-2.-3l^il^) 



m 4m 3m 



{viV2V3\i'2) 



— ai — :!!L — al 

m 2m m 



{2 2 1 

2 (iTg -V2 + 2^?3 2^2 - ^^gl^/'i) + (V3 - V2 ^^3l^2) 

(4) 

where the ^^(e) and ts{e*) are symmetrized two-body t-matrices with the two-body subsystem 
M _ M and e* = i? - |i - ^. 

4m dm 2m m 



energies e = E — ^ — ^ and e* = E — ^ — ^. The matrix elements of the 3BF term, 



U1U2 u 3\W{23\'^) , are evaluated in the next section. 



III. THE EVALUATION OF 3BF MATRIX ELEMENTS IN A 3D APPROACH 

Each part of a 3BF with two scalar meson exchanges and a constant meson-nucleon 
amplitude, which is shown in Fig. [H can be written in the following form 

^123«g2 + ^2 g/2+^2 

with a cutoff function 
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and momentum transfers Q and Q' 



Q 



ki - k 



1 




(123,4;12) 



{U2 - 'U2}(231,4;23) 



Q' 



k'2 - k2 




(123,4;12) 



{u'2 - M2}{312,4;31) 



(7) 



where the multiple indices for each curly bracket denote the two-body followed by the 3 + 1 
fragmentation. 

For the evaluation of Eq. (jlj) matrix elements of the form {u 1 

M2M3|W^iSl^) need to be 

calculated. From Fig. [1] we see that ^1^23 

can be considered as a sequence of meson exchanges 
in the subsystem (23), where it is called for convenience subsystem 1, and subsystem (31), is 
called 2. Since the structure of the 3BF we consider is specified by two momentum transfers 
of consecutive meson exchanges, it is convenient to insert a complete set of states of the 

type 2 between VF123 ^^d |\E') and another complete set of states of type 1 between the two 

(3) 

meson exchanges. Then the matrix element of W^i23 is rewritten as 



Here the subscripts 1, 2, 3 of the bra and ket vectors and in integrals stand for the different 
types of three-body coordinate systems of (3 + l)-type fragmentation {ijk,4:;ij). Both 
meson-exchange propagators in the 3BF term only depend on the momentum transfer in a 



3{UiU2U3\Wl2m 




X2{u'ru'^'u'^'\'^) 



(8) 
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two-body subsystem, as indicated in Eq. ([7]), i.e. 

I ->// -'II ->ll\ 



Q"^ + 171"^ 



jr6/-'i -'ii\i:it-'i -'ii\ 

5 [u^-u^)6 (M3 - u^)—j 



II, F{{u',-u'J,f) 



u'2 — U2Y + ml 



/ 2^ 



'Q'2 + ^2 



Fl ^2 ^3 /2 



'u'i" -u'''Y + m'i 



XO/^lll ^ llll\ / -> III -^llll\ 



Using Eq. ([9j), one can rewrite Eq. ([8]) as: 



(9) 



3(MiM2M3|W^12ll^) 

D^m' 2,{UiU2U'i\u\u'2u'^)l 



(n'2 - w'a')^ + ml 



X, / 7~l3 III i->i -fii ->i I ->lll ->lll ->lll\ 



3„,/w 



I -'III -'I 11 ^"'liTtV 

X 2\Mi ^2 "^3 1^/ 



(10) 



We would like to point out that in our vector based method the calculation of the trans- 
formations from one three-body subsystem to another, i.e. 3( | )i and i( | )2, are efficiently 
five-dimensional interpolations, whereas in calculation of the coordinate transformations via 
a PW decomposition, there is a complicated angular momentum recoupling algebra involved. 

Also we would like to mention that we do not follow the explicit evaluation of the co- 
ordinate transformations in Eq. ffTOl) leading to expressions with meson propagators which 
contain linear combinations of three or four momentum vectors. Thus direct integrations for 
evaluating the matrix element of the 3BF would involve magnitudes of momentum vectors 
and angles between all of them, which can be very complicated and involved. We therefore 
follow the method proposed in Ref. 33| and do not carry out the coordinate transformation 
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analytically, we evaluate the integration of Eq. (fTOj) in separate steps where in each step we 
only integrate over one vector variable at a time. Thus we define an auxiliary function 



771 / -> /// ^ /// ^ /// 



nil 

2 



U 



[u: 



u 



+ mi 



X 



1^) 



(11) 



the integration of the meson exchange between particles 2 and 3 in Eq. ffTTj) is carried 
out completely in the coordinate system of type 2. Once F2{u'l' ,u'^') is obtained, it 
needs to be expressed in terms of momenta in a coordinate system of type 1 in order to 
carry out the integration over the remaining meson exchange. This transformation, labeled 
F2i{u[,U2,u'^) is explicitly given as 

F2i{u[,U2,u'^) 

7-l3 /// /->! -> II ->l \ -> III -•III ->lll\ 17 ->/// ->lll\ 



'2 

F2i 



-u 



4 



-^11 ->i 

U2,Ul 



1 



U2,U'^) 



(12) 



Here we used that F2{u'i ,U2 ,11'^') is a scalar function due to the total wave function 
\E'('Ui -^2 %) being a scalar in the ground state. In our vector based method, this trans- 
formation is effectively a five dimensional interpolation on F2 in Eq. ffTTl) . which can be 
handled by the cubic Hermitian splines of Ref . 4l| . The integration over the second meson 
exchange between particle 3 and 1 in the coordinate system of type 1 is now given by 



Fi{u[,U2,u's] 



niu'2 



u 



ll\2\ 



il\2 



+ m; 



:}i F2i{u[,U2,u'3[ 



(13) 

'2-«2j^ + "^r • ^ - - 

The matrix element ^{ui U2 "^31^^1231^) is finally obtained by integrating Fi{u[, U2, 'U3) over 
u[, u'2 and Mg, i.e. carrying out the final coordinate transformation from the system of type 
1 back to the one of type 3, 

D^u' 3(wi -U2 Mal-u'i -u'a -Ug)! Fi(-u'^, ■U2, ^3) 

(14) 



Fi[--Ui - -U2, Ui - -U2, Us) 



9 



IV. CHOOSING THE COORDINATE SYSTEMS 



In order to solve the coupled three dimensional Yakubovsky integral equations, Eq. (jlj), di- 
rectly without employing PW projection, we have to define suitable coordinate systems. The 
Yakubovsky components are given as a function of Jacobi momenta vectors and as a solution 
of integral equations. Since we ignore spin and isospin dependencies, the both Yakubovsky 
components are scalars and thus only depend on the magnitudes of Jacobi momenta and 
the angles between them. The first important step for an explicit calculation is the selection 
of independent variables. As indicated in Ref. |331] one needs six variables to uniquely spec- 
ify the geometry of the three vectors. The coupled three dimensional Yakubovsky integral 
equations, Eq. (jl]), with only 2BFs was solved successfully in Ref. 40|]. For the evaluation 
of the 3BF term in the first Yakubovsky component in Eq. (jl]), z{uiU2 u ■i\Wi22,\^) we start 
with calculating first F2{u'i ,u'2 ,u'^'), Eq. (fTTj) . and realize that for this integration we can 
choose Ug" parallel to the 2;-axis and u'2 in the x — z plane. This leads to the simplification 
of the azimuthal angles. The explicit expression is 



r2{Ui , U2 , , , X2 , X^m^m) 

1 j'2lT 

dx'^' / #'2'" r«',<,y2"'2"") 

Jo 

y:\iJ di'" 11"" u'" t'" t"" r"^" 

^ ^ V""! 1 ""2 1 ""3 ) "''I ) -^2 ' •^u'-l'u2 



(15) 



with 





= lu'l'l 




III 

U2 


1 -'lll\ 

= F2 1 




III 

U3 


1 -'lll\ 

= Fa 1 




III! 

U2 


= 




X 


'-III '-III 

= % .u^ = 


cos(^9;") 


X2 


'-III '-III 

= % ■U2 = 




X2 




cos(??r 
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Hi'" 2"" = u'l ■U'2 



^ /// ^ //// 

y2"'2"" = ^2 •'U2 



+ ^1 - x-y 1 - xr cos(^r ) 



3 

J'l^,'" 



yii//2'» — x'"x2 



v/i-<Vi-4"' 

2/i"'2"" - x'/'za" 



3 



"1 "2 /1 /1 ^;w72 



2 



= «2 + ^2 + ^"^2 "^2 y2"'2"" 
r(„"" „ \ - -^(^r><^l/2'"2"") ..px 

J- 1^2 '"2 ■,y2"'2"") ^^fT- — 5 

Similarly for the integration over the second meson exchange, i.e., the evaluation of 
Fi{u[, u'2, u'^) of Eq. (1131) . we can choose u'^ parallel to the z-axis and m 2 in the x — z plane. 
This leads to the explicit expression which is functionally the same as Eq. (fT5|) : 

rx\U^, ^2, %, Xy, X2, ^u'^u'J 

f + l f2-K 

II 11 2 

VbV ■ 





rfWsVa / dx-2 \ #2 r(M2,M2,2/2'2"^ 

J -\ Jo 



2l(w'i, M2, W3, x[, X2, ^J^u'^) (17) 

with the same variables as Eq. f|T6|) with w'^, Wg, M3, Ug, x'^, X2, X2, (/^'i, V52 instead of 
<',<,<,<',x'/',x'2",xf (^f . The evaluation of F2i{u[,u'i,u'^), Eq. ([I2D, is not an 
integration but rather a five dimensional interpolation and explicitly is given by 

-^21(^1, "^2, %, X^, X2, X^?^,) 

= F2(ni,n2,M3,xni„^,xn2„^,Xn'n2) (18) 
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with 



1 3 1/9 
Hi = \ - -u[ - -u'^\ = -duf + -wf + 3m;<2/i'2" 



Ha = \u[ - ^u'^l = \ju'l + ^uf - u[u'^yr2" 

1 ^ ^3 1 1 3 
/ ->/ -> //\ ^ ' / II ^ II ii\ 

XUiu'^ — {—2'^l — ^U2)-U^——{ — -UiX^ — -U2X2) 
I -> I ^ -> ii\ " I ^ / I I ^ II li\ 



niHa' 2 ' 8 ^ 2 
^niHa - n——2 — — 



ni«3V 

Finally, the matrix element 3(m 1 ^2 'tis 1 1^123!^) is explicitly obtained by a five dimensional 
interpolation as 

(20) 



with 



, 1^ 3^ , 1/2 9 2 ^ 

^6 = I - ~ 4^2! = 2 Y "1 + + 3Miti22/l2 



1 ^ , / n 1 



TTy = \Ui - -M2I = y Mf + -M2 - UlU2yi2 

. iZ "3^ ^1^1 3 ^ 

/ 4 VTg Z 4 

XnrU3 = [Ul - 7;^2>-U3 = —{UiXi - -U2X2) 

2, TX-j 2, 

X^^^j = {--Ul - -U2).{Ui - -U2) 
1.12,32 1 



o^l + 0^2 - 7T"1^2?/I2j 



TTeTTr' 2 ' 8 " 2 

The last term of first Yakubovsky component in Eq. (jlj) requires an additional integration of 
the matrix element [ui U2U3\w'^^\^) and the fully off-shell two-body t- matrix. Again, with 



12 



choosing Ms parallel to the 2;-axis we only have four vectors to consider, Ui,U2,U3 and u'l, 
thus the integration is of a similar type as the one of the first three terms of first Yakubovsky 
component in Eq. (jl]), 



with 



^ Am 3m 



^ -'0 ^-1 ^0 E 1 T^ — ^ 

m 4m 3m 

(3), 



x({i;U2n3|M^123l*) (22) 



yyi = u[.ui = x'lXi + \J 1 — x'l \J 1 — x\ cos{lp'i — Lpi) 



yi'2 — u[.U2 = x[x2 + y 1 — xfy 1 — xl cos((/9'J (23) 

These considerations lead to the explicit representation for the Yakubovsky components 
and |?/'2): 
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U2 % Xi X2 X 



1 



m 4m 3m 



oo /■+! /.27r 



du2U2 / / (iy92 ^^^^(wi, TT, a;; e) 

J~i Jo 



+z/'l(7ri Tl2 VTs X22 2:23 X 



TVS 



+1p2{T^l 71-4 TTg X32 X33 X^ ) 



'01 



m 4m 3m 



xFi{nQ, tTy^u^, Xj^'^u^, ^J^^'^) 



1 

V^2Kf2t^3^1^2^12) = 72 

E — — i!2_ _ ]h 
m 2m m 

dv'^v'n dX'J d(l)',tsiv,,v',,Y,,r,e*) 

J-1 Jo 



X I 2^1(^73 Si S2 X12 Xi3 XjsJ 

+MV3 V2 v', X22 X23 X^'lJ I (24) 



The coupled equations, Eq. ( 12^ . is the starting point for numerical calculations and the 
details will be described in the next section. The 3D representation of total wave function 
1^) which directly appears in Eqs. f|T5l) and (124|) is represented in Ref. 40|, where we have 
presented it as function of vector Jacobi momenta. 

In a standard PW representation Eq. (jl]) is replaced by two coupled sets of a finite 



number of coupled integral equations 20], where the evaluation of two-body matrices and 
permutation operators P, P and P34 as well as coordinate transformations due to considering 
angular momentum quantum numbers instead of angle variables leads to more complicated 
expressions in comparison to our 3D representation. 
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V. NUMERICAL TECHNIQUES 



In this section we describe some details of the numerical algorithm for solving the coupled 



Yakubovsky three dimensional integral equations, and more details can be found in Ref. |40|. 
The Yakubovsky components are given as a function of Jacobi momenta vectors and as a 
solution of coupled three dimensional integral equations, Eq. (124|) . The both Yakubovsky 
components are scalars and thus only depend on the magnitudes of Jacobi momenta and the 
angles between them. The dependence on the continuous momentum and angle variables 
should be replaced in the numerical treatment by a dependence on certain discrete values. 
For this purpose we use the Gaussian quadrature grid points. The coupled Yakubovsky 
equations represent a set of three dimensional homogenous integral equations, which after 
discreatization turns into a huge matrix eigenvalue equation. The huge matrix eigenvalue 
equation requires an iterative solution method. We use a Lanczos-like scheme that is proved 



to be very efficient for nuclear few-body problems [42|. The momentum variables have to 
cover the interval [0,oo]. In practice we limit the intervals to suitable cut-offs and their 
values are chosen large enough to achieve cut-off independence. The functional behavior 
of the kernel of eigenvalue equation is determined by the two-body t— matrices. We also 
solve the Lippman-Schwinger equation for the fuUy-off-shell two-body t— matrices directly 
as function of the Jacobi vector variables [2^. Since the coupled integral equations require 
a very large number of interpolations, we use the cubic Hermitian splines of Ref. j41| for 
its accuracy and high computational speed. It should be mentioned that by adding the 
additional grid points, to all momentum and ±1 to all angle grid points, we avoid the 
extrapolation outside the Gaussian grids. 

VI. NUMERICAL RESULTS 

A. Three- and Four-Body Binding Energies 



In our calculations for 2BF we employ the spin-averaged Malfiiet-Tjon V potential 43|. 



This force is a superposition of a short-ranged repulsive and long-ranged attractive Yukawa 



interactions. We use the same parameters as given in Ref. 40| where the nucleon mass 
is defined by ^ = 41.470 MeVfm^. With this interaction we solve the Lippman-Schwinger 
equation for the fuUy-off-shell two-body t— matrices directly as function of the Jacobi vector 
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TABLE I: Three-body binding energies with and without three-body forces in MeV. The numbers 
in parenthesis are binding energies calculated in Ref. for three-body bound state with a modified 
version of Malfliet-Tjon by a cutoff function of dipole type. Also the number in bracket is calculated 
in FY scheme in PW representation, Ref. [l^. 



Potential 


Three-body Binding Energy 


MT-V 


-7.74 [-7.73] 


MT-V+MT3-I 


-8.92 


MT-V+MT3-II 


-8.70 


MT2-II 


-7.69 (-7.70) 


MT2-II+MT3-I 


-8.87 (-8.87) 


MT2-II+MT3-II 


-8.64 (-8.65) 



TABLE II: Four-body binding energies with and without three-body forces in MeV. The number 
in bracket is binding energy calculated in FY scheme in PW representation, Ref. ^19i] . 

Potential Four-body Binding Energy 



MT-V -31.3 [-31.36] 

MT-V+MT3-I -38.8 
MT-V+MT3-II -37.5 



variables as described in Ref. 28| . The so obtained t— matrices are then symmetrized to get 
ts{ui,7r,x;e) and ts{vi,v'^,Y^^'] e*). 

For four-body (three-body) binding energy calculations thirty (forty) grid points for Ja- 
cob! momentum variables and twenty (thirty two) grid points for angle variables have been 
used respectively. As demonstrated in tables [T] and [Tll the calculations of the three- and 
four-body bindin g en ergies using only the MT-V 2BF yield the values E = —7.74 and 



-31.3 MeV, Ref. 



40]. 



In our calculations for 3BF we use a model of 3BF which is based on multi-meson ex- 
changes. We study two different types of 3BFs, a purely attractive and a superposition of 



attractive and repulsive, which are named MT3-I and MT3-II respectively Ref. 33|. As 



shown in Ref. 



331 ] The parameters of these 3BFs are chosen so that the correction due 



to these 3BFs to the three-body binding energy calculated with the modified Malfiiet-Tjon 
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2BF (MT2-II) is small, and they lead to binding energies near to the triton binding energy. 

The three- and four-body binding energies calculated in 3D approach are given in tables 
HI and ini Our results for three-body binding energies with the addition of the MT3-I and 
MT3-II 3BFs, while MT-V is used as 2BF, are -8.92 and -8.70 [MeV] and while MT2-II is 
used as 2BF are —8.87 and —8.64 respectively. Our results agree with corresponding values 
presented in Ref. [33| and [29|. Our results for four-body binding energies with the addition 
of the MT3-I and MT3-II 3BFs, while MT-V is used as 2BF, are -38.8 and -37.5 [MeV] 
respectively. Unfortunately we could not compare our results for four-body binding energies 
with other calculations, since to the best of our knowledge no comparable work with scalar 
two-meson exchange 3BFs exists. So in order to test the accuracy of our calculations we 
carried out two numerical tests which are presented in next section. 

According to our experience for four-body bound state calculations with 2BF alone 40| , 



we expect that our results with 3BF provide the same accuracy in comparison to other cal- 
culations of the four-body binding energy based on PW decomposition, while the numerical 
procedure are actually easier to implement. 



B. Test of Calculations 



In this section we investigate the numerical stability of our algorithm and our 3D repre- 
sentation of Yakubovsky components. We specially investigate the stability of the eigenvalue 
of the Yakubovsky kernel with respect to the number of grid points for Jacobi momenta, 
polar and azimuthal angle variables. We also investigate the quality of our representation of 
the Yakubovsky components and consequently wave function by calculating the expectation 
value of the Hamiltonian operator. 

In table Hm we present the obtained eigenvalue results for binding energies given in tables 
[T] and [Til for different grids. We choose the number of grid points for Jacobi momenta as 
Njac, for spherical angles as Ngph and for polar angles as Npoi. As demonstrated in this 
table, the calculation of the eigenvalues A convergence to the value one for Njac = 30 and 
Ngph = Npoi = 20. It should be clear that the solution of coupled Yakubovsky three- 
dimensional integral equations, with six independent variables for the amplitudes, is much 
more time-consuming with respect to the solution of three-dimensional Faddeev integral 

n 

equation [33], with three variables for the amplitude. 
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TABLE III: Stability of the eigenvalue A of Yakubovsky kernel with respect to the number of grid 
points in Jacobi momenta Njac-: spherical angles Nsph and polar angles Np^i- Emt-v = —31.3, 
Emt-v+mtz-i = —38.8, Emt-v+mt3-ii = —37.5 MeV and Ai,A2 and A3 are corresponding 
eigenvalues. 



Njac 


Nsph = Npoi 


Ai A2 A3 


20 


20 


0.987 0.988 1.010 


26 


20 


0.995 0.996 1.004 


30 


12 


0.997 0.997 1.003 


30 


16 


0.999 0.999 1.001 


30 


20 


1.000 1.000 1.000 



The solution of coupled Yakubovsky three-dimensional integral equations in momentum 
space allows to estimate numerical errors reliably. With the binding energy E and the 
Yakubovsky components l-^i) and 1-02) available, we are able to calculate the total wave 
function from Eq. ([2]) by considering the choice of coordinate systems which are rep- 



resented in Ref. 40|. So in order to demonstrate the reliability of our calculations we can 
evaluate the expectation value of the Hamiltonian operator H and compare this value to 
the previously calculated binding energy of the eigenvalue equation, Eq. fl24|) . Explicitly we 
evaluate the following expression: 

= 12 {i;,\Ho\^) + 6 {^P2\Ho\^) 

+ 4(^|W^i23|*) (25) 

where V represents the 2BFs Vij) and W the 3BFs iJ2i<j<k ^ijk)- The expectation 

value of the kinetic energy (Hq) and the 2B potential energy (V12) have been evaluated in 
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Ref. |40j. The expectation value of the 3B potential energy, (VF123), is given by 



1 1^1231^) 




3 X Stt^ / dui u\ / dxi i difi 



Jo J -I Jo 




X 



/ du2 ul / dx2 / dus ul 
Jo J~i Jo 



X 



\E'('Ui U2 U3 Xi X2 V^l)W^i23 U2 U3 Xi X2 V?l) 



(26) 



Here the integrations need the evaluation of the matrix element 3(n 1 u 2 u 3|iy{23|^) of 
Eq. ( JTOl) . The expectation values of the kinetic energy (Hq), the 2B interaction {V), the 3B 
interaction (W) and the Hamiltonian operator (H) for three- and four-body bound states 
are given in tables IIVI and |V] respectively. In the same tables the corresponding binding 
energies calculated in 3D scheme are also shown for comparison to the expectation values of 
the Hamiltonian operator. One can see that the energy expectation value and eigenvalues 
E agree with high accuracy. All these numbers are not meant to provide insight into the 
physics of three and four interacting nucleons, but serve only as a demonstration that this 
technique allows a very accurate and easy handling of typical nuclear forces consisting of 
attractive and repulsive (short range) parts. In addition, they will serve as benchmarks for 
future studies. 

3BF effects have a stronger impact on four-body bound state than for three-body bound 
state as can be seen for instance by comparing expectation values of the potential energies 
for the two systems. We find in case of MT-V (V) = —37.51 (—101.0) MeV for three 
(four)-body bound state without 3BF and {V) = -40.63,-40.02 (-110.1,-107.6) MeV 
with MT3-I and MT3-II 3BFs correspondingly. In the latter case the expectation values 
for the 3BFs are (W) = -1.41, -1.07 (-7.5, -6.0) MeV for three (four)-body bound state. 
Already the trivial fact that there are four triplets in four-body bound state makes it clear 
that one has to expect 3BF effects to be more pronounced in the four-body bound state 
than in the three-body bound state. 
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TABLE IV: Expectation values with respect to the three-body wave functions for various potential 
combinations. We present the expectation values of the kinetic energy (Hq), the 2B interaction 
(V) and the three-body interaction {W). Additionally the expectation values of the Hamiltonian 
operator (H) are compared to the binding energy results from the Faddeev equations. All energies 
are given in MeV. 



Potential 




(V) 


{W) 


(H) 


E 


MT-V 


29.77 


-37.51 




-7.74 


-7.74 


MT-V-FMT3-I 


33.13 


-40.63 


-1.41 


-8.91 


-8.92 


MT-V-FMT3-II 


32.38 


-40.02 


-1.07 


-8.71 


-8.70 


MT2-II 


28.64 


-36.33 




-7.69 


-7.69 


MT2-II4-MT3-I 


31.88 


-39.40 


-1.34 


-8.86 


-8.87 


MT2-II-FMT3-II 


31.17 


-38.78 


-1.04 


-8.65 


-8.64 



TABLE V: The same as table IIVI but for four-body case. 
Potential {Ho) {V) {W) {H) E 
MT-V 69.7 -101.0 - -31.3 -31.3 
MT-V-FMT3-I 78.8 -110.1 -7.5 -38.8 -38.8 
MT-V-hMT3-II 76.1 -107.6 -6.0 -37.5 -37.5 



VII. SUMMARY AND OUTLOOK 



Instead of solving the coupled Faddeev- Yakubovsky equations in a partial wave basis, we 
introduce an alternative approach for four-body bound state calculations which implement 
directly momentum vector variables. We formulated the coupled Yakubovsky equations for 
identical spinless particles, interacting by two- and three-body forces, as function of vector 
Jacobi momenta, specifically the magnitudes of the momenta and the angles between them. 
We expect that coupled three-dimensional Yakubovsky equations for a bound state can be 
handled in a straightforward and numerically reliable fashion. In comparison to an angular 
momentum decomposition which is commonly used [l9|-|23], this direct approach has great 
advantages. In our Three-Dimensional case there is only two coupled three-dimensional in- 
tegral equations to be solved, whereas in the partial wave case one has two coupled sets of 
a finite number of coupled equations with kernels containing relatively complicated geomet- 
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rical expressions. The comparison of 3D and PW formalisms shows that our 3D formahsm 
avoids the very involved angular momentum algebra occurring for the permutations and 
transformations and it is more efficient especially for the three-body forces. 

The three dimensional Yakubovsky integral equations was successfully solved using 
Malfliet-Tjon type 2BF alone, and its numerical feasibility and accuracy established 40|. 
Here we present results including the scalar two-meson exchange three-body force and study 
its effects on the energy eigenvalue and the four-body wave function. The stability of our 
algorithm and our Three-Dimensional representation of Yakubovsky components have been 
achieved with the calculation of the eigenvalue of Yakubovsky kernel, where different num- 
ber of grid pints for Jacobi momenta and angle variables have been used. Also we have 
calculated the expectation value of the Hamiltonian operator. This test of calculation repre- 
sents good agreement between the obtained eigenvalue energy and expectation value of the 
Hamiltonian operator. 

This is very promising and nourishes our hope that calculations with realistic two and 
three-nucleon forces, namely considering spin and isospin degrees of freedom, will most likely 
be more easily implemented than the traditional partial wave based method. 

To this aim the first step for realistic calculations of three-nucleon bound state in a 



realistic Three-Dimensional approach 



ing energy with Bonn-B potential 



las already been taken by calculation of Triton bind- 



44 



45l ] and formulation of four-nucleon bound state is 



currently underway and it will be reported elsewhere 46|. They will be the first steps for 
realistic calculations of three- and four-nucleon bound states in a Three-Dimensional scheme. 

It should be mentioned that the input to such calculations is the NN t-matrix which is 
calculated in an approach based on a helicity representation and depends on the magnitudes 



of the initial and final momenta and the angle between them [31|- Consequently the cal- 
culation of NN t-matrix in helicity representation needs the NN potentials in an operator 
form which can be incorporated in 3D formalism. As indicated in sec. 3.2 of Ref. [47| (or 
sec. HI of Ref. [3l|) the general structure of the NN potential operator which fits well to 



the helicity representation is given, and on this representation both Bonn-B and AV 



potentials are given in operator form, see appendixes C and D (or sec. IV of Ref. 3l|). 



8 NN 
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